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t Abstract: I propose to consider photon tunneling as a space-time correla- 

■ tion phenomenon between the emission and absorption of a photon on the two 
| sides of a barrier. Standard technics based on an appropriate counting rate 

formula may then be applied to derive the tunneling time distribution without 
any ad hoc definition of this quantity. General formulae are worked out for a 

■ potential model using Wigner-Weisskopf method. For a homogeneous square 
barrier in the limit of zero tunneling probability a vanishing tunneling time is 

• obtained. 
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§ ■ 1 Introduction 

The phenomenon of tunneling was recognized immediately after the birth of 
quantum theory as one of the most striking features of microphysics. Though 
i tunneling probabilities are predicted unambiguously by the theory this is not 

<*S || the case with the time required for tunneling 1 . Tunneling time is equal to the 

+Jj ■ time interval required to travel a certain distance when a barrier between the 

Ch I endpoints is present minus the time required to travel the distance outside the 

barrier. This quantity is directly measured in an experiment "with clocks" i.e. 
qh using two detectors that signalize the moment of departure and the moment of 

arrival of the tunneling particle (see Fig.l), but no experiment of this type has 
. so far been actually performed. No theory is required to infer the value of the 

. tunneling time from such an experiment because it is just this experiment which 

defines the notion " tunneling time" operationally. It may of course happen that 
an experiment if performed would give no sharp tunneling time but rather a 
tunneling time distribution. 

Theory is needed to predict the result of the experiment. However, up to now 
no method has been proposed which would permit us to infer from quantum 
theory the expected value (or distribution) of a tunneling time experiment as 
described above. The Wigner-time 0, the Biittiker-Landauer time [3], 0] and 
the Larmor-time [5] are each based on different secondary criteria for the time 
spent by the particle under the barrier but neither of them has been shown to 
be equivalent to the primary notion of the tunneling time operationally defined 
by an experiment with clocks. The problem has recently been surveyed by 
R.Y.Chiao (see also [7], 0). 

1 E-mail: peter@hrasko.com 
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The experimental investigation of the problem became possible in the last 
decade thanks to the recognition that massive particles may be replaced by 
photons. In the experiments with photons periodic layers of alternately high 
and low index media serve as optical tunnel barrier. Such experiments, while of 
great importance in themselves, may shed light on the tunneling of nonrelativis- 
tic massive particles as well, because the mathematical form of the equations, 
governing the tunneling process, are essentially the same for both cases. 

The great advantage of using photons consists in the possibility of converting 
time intervals into phase shifts. This trick, going back to the Michelson-Morley 
experiment, makes the measurement of exceedingly small transit times feasible 
by the shift of the interference pattern in an appropriate interferometer. In 
the Berkeley-experiments [5], ^21 UV photons were split into a pair of photons 
of equal frequency by spontaneous frequency down conversion in a nonlinear 
medium. In a Hong-Ou-Mandel spectrometer the two photons produce charac- 
teristic minimum in their rate of coincidence when pass through a beam-splitter 
which is in symmetrical position with respect to the beams. When, however, 
a thin tunneling layer was posed into the path of one of the beams the coinci- 
dence minimum could be maintained only if the beam-splitter was displaced at 
a certain distance. From this distance the tunneling time of the photon could 
be deduced and turned out to be about 2 fs while the time required to traverse 
the same distance in vacuum is 3.6 fs. 

This result which is in rough agreement with the Wigncr tunneling time 
indicates that the speed of tunneling may exceed the speed of light in vacuo 
by about 70 percent. It may be noted, however, that stationary interference 
experiments, contrary to the experiments with clocks, permit us to infer tunnel- 
ing time only at the expense of imagining photons, moving in the arms of the 
equipment, as particles. It is this picture which might suggest the counterfac- 
tual conditional statement that if tunneling time was directly measured it would 
be equal to the time inferred from the stationary interference experiment. In 
quantum physics, however, counterfactual statements based on imagined rather 
then real events are in general invalid even if their validity in classical physics 
would be beyond doubt. 

This difficulty has been clearly recognized in jS], where it was pointed out, 
that different experiments may lead to different tunneling times. Based on this 
fact the suggestion was made to abandon the unique definition of the tunneling 
time and replace it by a multitude of equivalent notions, assigning to each 
conceivable experimental setup its own value of this quantity. In particular, 
tunneling times measured in an experiment with clocks and deduced from a 
particular stationary experiment are both legitimately called tunneling time 
albeit under different circumstances. 

But relativization of concepts often leads to confusion. It seems, therefore, 
safer to keep consistently to the definite meaning of tunneling time as formu- 
lated at the beginning of this section, because it conforms with the general use 
of the word " time" , including its operational meaning. In doing so one naturally 
has to admit, that the time parameters inferred from experiments of other types 
are in general different from what is properly called tunneling time. This point 
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of view may, perhaps, be argued for by noticing that in relativity theory only 
such velocities are of significance which may be related to real pairs of events 
through the elementary formula As/ At while velocities inferred from station- 
ary experiments are necessarily based on imagined rather than real events. It 
may indeed be dangerously misleading to call " velocity" a quantity which lacks 
essential connotations of this term. 

The purpose of the present work is to suggest a theoretical scheme to cal- 
culate the result of an experiment with clocks. The main obstacle on the way 
to formulate such a method is the lack of a rigorous quantum theory of the 
arrival time distributions of photons and particles. In quantum optics this diffi- 
culty has been overcome by the replacement of the statement "The detector has 
clicked" with the statement "The detector atom is in one of its excited states". 
Based on the assumed equivalence of these two statements (see Section 8) work- 
ing formulas, known as counting rate formulas, were derived to first order in 
the photon-detector interaction which have since been used succesfully in 
treating space-time correlations between photons. 

The direct measurement of the tunneling time of photons is actually the 
measurement of the time correlation between the emission and the absorption 
of the photon on different sides of the barrier. If, therefore, one accepts the basic 
assumptions underlying counting rate formulas, an appropriate formula of this 
kind can be derived which permits us to calculate the value (or distribution) of 
the tunneling time as operationally defined. 

In the present work we confine ourselves (1) to the presentation of this 
formula and (2) to show that in the absence of any barrier it leads to the expected 
time correlation. In addition, it will be illustrated on a simplified model, how 
the theory works when a barrier is present. No attempt is made to apply the 
method to realistic barriers as e.g. to that used in the Berkeley-experiment. 

For a very high and broad barrier our model calculation gives sharp tunneling 
time which is equal to zero and so it agrees qualitatively with the result of the 
Berkeley-experiment which predicts the reduction of the time required to travel 
a given distance when a barrier is present. This preliminary result may be an 
indication that the time parameter inferred from the Berkeley-experiment is 
indeed closely related to the tunneling time as defined in an experiment with 
clocks. 

2 Working formula for a thought experiment 
with clocks 

The experimental setup is shown on Fig. 1. The photon source will be a two- 
level atom at the origin of the coordinate system. An atomic detector at the 
point r = (0, 0, z) of the z coordinate axis serves to detect the photon. The 
barrier, an infinite homogeneous layer, supporting evanscent waves in a broad 
interval around the photon wavelength, is placed between the source atom and 
the detector perpendicular to the z-axis within the region (a, b) of this axis. An 
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appropriate source detector is assumed to be present in the immediate vicinity 
of the source atom which clicks at the moment of the photon emission. It may 
be assumed that it detects a particle (or another photon) which accompanies 
instantaneously the emission of the tunneling photon. In short, the setup works 
on the principle of the usual time of flight spectrometers. 

Delay time analyser 



Source detector 



Source O- - 



Photon 




Optical 
barrier 



D 



Photon 
detector 



Figure 1. 



Assume that at t = the source atom and the detector atom are in their 
excited and ground states respectively and no photon is present. Our aim is to 
calculate the probability density w(t\,t<i) of the source detector clicking at t\ 
and photon detector clicking at t 2 . 

In order to calculate w(t\,t2) one has first to determine the probability 
p(t\,t2) of finding the source in its ground state at the moment t\ and the 
photon detector in some of its excited states at ti- Both t\ and t% are chosen 
arbitrarily "by ourselves" rather than by the experimental setup itself. Then, 
having p(t\, t 2 ) calculated, w(t\,t2) is obtained by differentiation: 

. x d 2 p(t 1 ,t 2 ) 

This kind of detour through p(t\ 1 t2) has been used since the sixties for the 
calculation of photon space-time correlations in light beams prepared in various 
ways. In first Born-approximation the detector degrees of freedom can be easily 
eliminated and p(t\, t 2 ) and its generalizations can be expressed in terms of the 
expectation value of some product of field variables. Relations of this type are 
called counting rate formulae The restriction to first Born-approximation 
ensures the nonnegativity of w(t\,t2)- 

This same procedure will be adapted below for the calculation of the time 
correlation in photon tunneling. The version of the counting rate formulae 
appropriate for the purposes of the present work is the following: 

/•oc 

p(t 1 ,t 2 )= dwflr(w)-|M w (ti,t 2 )| 2 , (2) 
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where a(ui) is the spectral sensitivity of the photon detector and 

M u {t u t 2 )= ^ dt-e^t&,va^T{^{ttfP?(t{))\e,vac). (3) 
Jo 

The amplitude in the last integral is the matrix clement between the states 
\g,vac) = \g) <g> \vac) and \e,vac) = |e) £g> \vac) of the photon vacuum and the 
source in its ground and excited states. P g is the projector to the ground state 
of the source atom, tp is the field operator of the photon, assumed spinless, and 
r is the position of the photon detector. The superscript H indicates Heisenberg 
picture in which the dynamical variables are driven by the Hamiltonian H = 
H s + Hf + H s f = Hq + H s f which is the sum of the source, the field and the 
source-field interaction Hamiltonians. The last term is assumed to be of the 
simple form 

H sf = Q-i P ('r = 0). (4) 

The operator Q acts in the Hilbert space of the source, its nonzero matrix 
elements being Q eg = Q* ge - The symbol T means time ordering. The main 
steps of the derivation of and © are summarized in Appendix A. 



3 Quantization of the photon field 

The equation satisfied by the field operator ip is 

§-^ + •^ = (5) 

(the system of units c = h = 1 is adopted). The last term represents the barrier, 
confined to the region (a, b) of the z-axis between the origin and the position 
of the photon detector at z. The width D of the barrier is, therefore, equal to 
b — a. For the time being no homogeneity of the barrier in the z-direction will 
be required. 

The axial symmetry of the experimental setup suggests that quantization be 
performed in cylindrical coordinates: 

V>M = E E f" ^7#^ \^r,^)U^r\r)v.{^\z)^e- i ^+ 
i Jo \/2uj 

m=—oo s=± u v 

+ a + s K,^)[/;„K|r)<(^|z)e-^e^] , (6) 
where 

W = V^r+^l (7) 

and 



U m (u! r \r) = ^- — J m {u) r r). (8) 
In 



The functions v±(u) z \z) obey the equation 

d 2 v± 



dz 2 



+ V{z)v±=u 2 z v±, (9) 
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The index ± indicates the direction of the incoming wave: 



V ^ = { Te i^(z-a) ifz>5j ( 10 ) 

j Te - iuJz (z-b) ifz<a? 

V-{UJ Z \Z) = \ ■ I u\ / l\ (11) 

\ e -w,[z-b) + R , e tu x {z-b) ifz> 5 7 

where we have used the fact that the transmission coefficient is independent of 
the incoming direction (T" = T, see Appendix B). These functions obey the 
relation 



E 



du z ■ v s (u) z \z)v*(lj\z') = 2ir ■ 6(z — z'). 



Equation |J3J is of the form of a Schrodinger-equation for a particle moving along 
the z-axis ([El, E3)- Therefore, the transmission and reflection coefficients T 
and R are the same which are found in the textbooks on Quantum Mechanics. 
One has only to notice that, as the asymptotic forms l|l(J|l . (|11|) indicate, T and 
R belong to a barrier of the given shape shifted to the origin (i.e. to a = 0, 
b = D). 



4 The Wigner-Weisskopf approximation 

In the next section the matrix element in J2J will be calculated in Wigner- 
Weisskopf (WW) approximation ^Hl- This approximation scheme is based 
on two assumptions. 

The first assumption is the confinement of the electromagnetic interaction 
to the subspace spanned by the vectors \e,vac) end \g, 1 photon). The special 
form of the interaction and the relation J m (0) = S m0 allow us to consider 
m = photons only. We have, therefore, the state vector of the form 

/>oo 

\t) =c{t)\e,vac) + y2 <L}rdu z Al, rU)z (t)\g,UrU}zs), (12) 
in the interection picture, which obeys the Schrodinger-equation 

^ = H sf (t)\t), (H sf (t) = e iH ^H sf {0)e- iH ^) , (13) 

and the the initial condition \t = 0) = |e, vac). It follows then that the coeffi- 
cients in H12|) satisfy the equations 

id (t) = J2 durduzA^^ty^-^^va^Hsfmg^rUOs) (14) 
«=± Jo 

iAl^(t) = c{t)e l ^-^ t {g^ r uj z \H sS {Q)\e,vac). (15) 
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In these equations f2 is the excitation energy E e — E g of the source, and u> is 
given by J7J). 

The second assumption of the WW-method consists in the Ansatz c(t) = 
e -Tt/2 p rom jj^j we ^ ave 

= / dt' ■ e 1 ^ - ° + ir / 2 )* / ( ff ,a; r t I ;,*|ff, / (0)|c ) t;ac) = 



~V ^ ^ z|0) c-» + zr/2 ■ (16) 

Now the value of T could be computed from l|14fl but since it is irrelevant for 
the present work we will not pursue this line any further. 

We notice that the exponential Ansatz is invalid for very short times of 
the order of 1/fi. This Ansatz prescribes the nondecay amplitude (e,vac\t) as 

eT w hose time derivative at t — is equal to —T/2. However, according to 
<|13|) . this time derivative actually vanishes: 



<9(e, vac\t) 
dt 



= (e, vac\H s f(0)\e, vac) = 0. 

t=o 



Owing to this shortcoming of the WW-approximation our subsequent consider- 
ations must be confined to the inside region of the future light cone of the state 
preparation event at r = 0, t = 0. To meet this condition the relations ii, t2 > z 
- or, more precisely, (ti — z)Cl ^> 1, (£2 — z)£l ^> 1 — will be assumed from 
the outset. This restriction is, however, irrelevant for a real experiment unless 
the latter is specially designed to investigate the outside region of the light cone 
(see Section 8). 

The interaction picture of the present section is connected with the Heisenberg- 
picture employed in @ by means of a unitary operator W(t,0): 

Q(t) = W(t, 0)O H (t)W+(t, 0) (17) 

(notice that the absence of an upper index on operators indicates interaction 
picture). Then, for (|12f> we have 

\t) = W(t,0)\e,vac). (18) 

Moreover, as a consequence of the first assumption the radiative corrections to 
\g,vac) must be neglected and we have the stability condition 

W(t, 0)\g, vac) = \g,vac). (19) 
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5 The formula for M^it^t?) in WW approxima- 
tion 

The formula © for the amplitude M u (ti, t 2 ) can be rewritten in the form 
M u> (t 1 ,t 2 ) = / dt-e lujt (g,vac\P"(t 1 )ip H (t,r)\e,vac)+ 



2 (2°) 
+9{t 2 -ti) / dt-e lujt (g,vac\[ V H (t,r),P g H (ti)]\e,vac), 
Jti 

where 9 is the step-function. 

When t\ > t 2 , the probability amplitude of finding the source deexcited and 
the photon detector excited is given by the first term of l|20() alone. One may 
expect that this amplitude must not depend on t\ since at the earlier moment 
t 2 when the photon detector was found excited the source was sure to having 
been already in its ground state. In WW-approximation, due to the stability 
condition (|19|) . this assertion is indeed true since t\ drops out of the integrand 
of the first term: 

{9,vac\Pg '{t{)ip H (t,r)\e,vac) = 
= (g 1 vac\W + {t ll O)P g (t 1 )W(t 1 ,0)W + (t,OMt,r)W(t,0)\e,vac) = 
= favaclPgitJWitutMtrflt) = (g,vac\^(t,r)\t). 

In addition to the stability condition use was made of the fact that, in the 
interaction picture, P g {t{) leaves \g,vac) invariant. We can, therefore, write 

MJt 2 ) = [~ dt-e iu,t (g,vac\P g H (t 1 )<p H (t,r)\e,vac) = 

( 21 ) 

dt-e lu}t (g,vac\p(t,f)\t). 

o 

Using (|12fl . the matrix element here can be cast into the form 

/>oo 

(g,vac\ip(t,r)\t) = V] / dw r d^zAl ;TU:s (t){g,vac\(p(t,T f )\g,^ r uj z s}. (22) 

s=± J 

The function A® u (t) is given by l|16|l while the matrix element in the integrand 
can be calculated using |JB}: 

(g,vac\f{t,f)\g,uj r UJ z s} = -^~Jy- e~ tujt v s (u) z \z) (23) 



ii 



(remember that r — (0, 0, z)). 
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The sum over s can now be performed: 

v*(oj z \0)v s (u z \z) = (e^ a + i?*e"^ a ) Te w ^ z ~ a ) + 
^ e -iu z (z - b) + R i e iu z (z - b)^ 

f rpft* e —2iuj z a _|_ rp* R / e —2iuj z b 



s=± 



According to (|48|l of the Appendix B the sum in the parentheses is equal to 
zero: 

TR * e -2iu z a + T * R > e -2w z b = 0; ( 24 ) 

therefore 

v* s (lo z \0)v s (lo z \z) = Te luJ * z + T*e~ luJ * z . (25) 

Putting now , ill'- it and i|25|l into (1221) and changing the integration variable 
from u> r to u> — WuJ^ + uj z we obtain 



(g,vac\P g (ti)ip (t,r)\e,vac) = (g,vac\tp(t,r)\t) 
did ■ ; / doj. 



(26) 

The matrix element in the second integral of l|2l)[l can be handled analogously: 
(g,vac\[<p H (t,?),lf (ti)]|e,«ac) = 
Y] / d^rd^z[Al rU)x (ti) - A s UrUz {t)]{g 1 vac\ip{t 1 r)\g,uj r uj z s) = 

s=± J ° 



Qge f -i(n-iT/2)t 



du) 



ii 



8?r 2 [ J u-tt + iT/2 

-»(n-tr/2)i! / ^ / d 



dw 2 



T(w z )e^ z +T*(^)e~^^ 
(27) 



Here t>t\ since this expression is the integrand in the second term of Q2U[). 
Comparing (|27|l with l|26|) we see that in WW-approximation 

(ff,mc|Pf (^"(i.rOMac) = ( ff ,«ac|[ V ff (t,f),i^(0)]|e,«ac). (28) 

6 Time correlation when no barrier is present 

In this case T = 1, the lo z integral in l|27|l gives 2sinajz/z and we have 

(g, vac\ [<p H (t, r), Pf (t x )] |e, mc)° = [if + 7° + 7° + 7°] , (29) 



where the superscript indicates the absence of the barrier and 



!0 iz) = ^ e -i(n-tT/2)t I duJ 



iz J uj — Q + iT/2 

1 vn wow f°° e -iw(t-z-h) 



iz J uj — £1 + iT/2 



l3 (z) = -ie-*(°-* r / 2 )* / <L> 



e -iuiz 



(30) 



iz 7 uj — Q + iT/2 

i2 W = i e -(n-T/2)tx / d w e 



Consider /f. The integration contour can be deformed upward to contain the 
positive imaginary axis and the quarter of the large circle at infinity. Since ui z 
is positive, the integrand is exponentially small on this part of the large circle 
and gives no contribution. Moreover, the integrand is regular in the upper half 
plane and no pole contributions arise. We have, therefore, 

This integral can be expanded in terms of the inverse of the large distance z, 
the leading term being of the order of 1/0,2, and so If turns out to be of second 
order. Therefore, in the leading (linear) order in 1/z, I® must be neglected. 

The same conclusion applies to the sum of 7° and 7° as well. In both of 
these terms pole contributions arise which are linear in 1/z but they drop out 
of the sum. 

Consider now I®. When t — z — t\ > the contour has to be deformed 
downward and a pole term 

/0 (z) = }_ e t t _ z _ ti) . 2m . e -i(n - $T/2)(t - z) (31) 
iz 

arises. In the opposite case of t — z — t\ < the deformation is upward and no 
pole is to be dealt with. The contribution of the integrals along the imaginary 
axis is negligible only if \t — z — ti|f2 ^> 1. Therefore, I®(z) is given by Ij31|l 
provided the step function is assumed smoothed on the scale 1/VL. From an 
observational point of view such a smoothing is of no significance and in what 
follows no attention will be payed to it. 
We have, therefore 

(g,vac\[<p H (t,r),P* (ti)]|e,mc)° = 

== «2£ ff(t _ tl _ z) e-i(n-a , /2)(*-»). (32) 

47TZ 
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From we obtain 

M° J (t 2 ) = - t dt- e iu \g, V ac\[y H {t^pX{Q)\\e,vac) = 



= j «9£ e(t 2 - z)e< n ~ lT / 2 ) z — . 

Attz y ' uj-tl + iT/2 

(33) 

The ^-function here may in fact be omitted since t 2 > z by assumption. 
For the integral, occurring in the second term of l|2(J|) we obtain 

„iu>t 



dt ■ e lujt (g, vac\ [y H {% r), Pf (0)] \e, vac) = 

i^t-t - */2)z e± - " + _ e ^ - n + l r/m 1 + z) 

47TZ — iZ + jI 1 2 



(34) 



Since 0(ti) = 0(t 2 — z) = 1, the right hand side of (|34[1 is equal to 

0(t a -ti-*)[M°(ti + z)-M°(f a )]. 

Hence, we have from (|20() 

M°(*i,* 2 ) = M° (ta) + 0(t 2 - ti - (*i + *) - M°M] = 

= 6{t x + z- t 2 )M° J (t 2 ) + 6{t 2 -h- z)M°(ti + z). 

Substituting this into 10 we obtain 

/>oo 

p°{t 1 ,t 2 )=6{t l +z-t 2 ) / duj-cr(u>)\M J {t 2 )\ 2 + 

Jo 

poo 

+0(t 2 -h-z) / du:-d{uj)\Ml{ti + z)\ 2 . 



(35) 



In arrival time measurements the spectral sensitivity must be as broad as pos- 
sible so we assume cr(a>) — a — constant. Then, substituting H33[) . we hnd 

jT 4, • fe)| 2 = • i (l - e"^ " *)) , (36) 

which leads to 

p o (tljta) = . i jx _ + z _ <2)e -r(t 2 - z) _ 0(t2 _ h z)e -Th I 

(37) 

The formula indicates that iu(ti, £2) may be different from zero only around 
t\ + z — t 2 = 0. Putting 

w°(t 1 ,t 2 ) = W{t 2 )-8{t 2 -t 1 -z), 
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wc find 



W(t 2 



t 2 -z+e 



dh ■w°(t ll t 2 ) = 



t 2 —z — e 

dp°(ti,t 2 ) 



dt 2 



tl=t 2 -Z + £ 



t 2 —z — t 

d P °(h,t 2 y 



dhdt 2 



dt 2 



ti=t 2 —z—e 



By (I37|) the second term is zero while in the limit of e = the first one is equal 
Airz 2 



Attz 2 



S(t 2 



z). 



(38) 



Though this is just the expected result it is far from being an obvious conse- 
quence of the counting rate formula @ , © an d of the reasoning in Appendix 
A which led to it . 



7 Time correlation in the presence of a rectan- 
gular barrier 

Let us choose in (0 a constant V(z) equal to fi in the interval (a, b) of width D 
and zero outside. Since is, up to constant coefficients, identical to the non- 
relativistic Schrodinger-equation, the transmission coefficient of this rectangular 
barrier can be taken over from quantum mechanics: 

■ n n l~2 2 1 - e Aia ^ UJ ^ 

T(u z ) = e- IU} * D -e- D VV 2 -u 2 z i — e (39) 

i _ Aia(w z ) . p -2D^fi 2 - uj 2 



in which 



UJ Z + iy/ [I 2 — UJ 



2 



The function T(uj z ) has a cut along oj z > fi. The physical values are those on 
the upper edge of the cut where y/ /j, 2 — lo 2 = —iy/uj 2 — /i 2 . 
The right hand side of l|27[) can now be written as the sum 

(g, vac\ [<p H (t, f), P g H (h)] \e, vac) = 9^[l 1+ J 3 + J 3 + J 4 ] , 
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where 

h{z) = -e-^-^l 2 )^ 
h{z) = e- i ( n - iV l 2 ) t 



lu -tt + iT/2 J Q 
e —iu(t — tt) 



du> 



LU-n + iT/2 Jo 
e —iw(t—ti) 



du z ■ T{u z )e iu) * z , 
<ko z -T*{Lo z y- luJ *\ 



du> 



u-fl + iT/2 



o 



dLu z -T*(Lu z )e- luJ * z . 

(40) 



We will argue below that the formula l|38[) . derived in the absence of a bar- 
rier, remains to reasonable accuracy valid also for a sufficiently high and broad 
rectangular barrier provided z is replaced in it by z — D. 

The variable z in the formulae l|4l)|l appears only in the inner integrals. If 



the first exponential factor of T is separated: T{u> z ) = e 



-iu> z D 



T(u) z ), then the 



exponentials in these integrals become e^^ii 2 

The essential contribution to the Ij-s must come from the region of integra- 

_on /.,2 _ 772 

tion around fl. If /i ^> Q, the exponent e t* M z i n the denominator is 

very small in this region. We expand the fraction in T(tu z ) in terms of this small 
quantity and retain from the resulting asymptotic expansion the first term only 
(i.e. we disregard the second term of the denominator). Then, 



1 - 



(41) 



Consider the behaviour of the factor e~ V" — * when \uj z \ — ► oo along 
some direction tp on the upper half o^-plane. Then the value ip — corresponds 
to the upper edge of the cut and approaching infinity we have yj fj? — ui 2 ~ 
— i\u) z \. Therefore, for positiv if 



so that 



-i\ui z \e ll P — — i\uj z \ ■ costp + \uj z \ ■ simp, 



approaches zero exponentially when \uj z 



Consider the inner integral 



F(u) 



dw z -T{u z )e iuj ^ z - D ) 



in 140(1 . Since its integrand vanishes exponentially on the large circle in the 
upper half plane the integration contour can be deformed in this direction into 
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two semiinfinite straight lines along the positive imaginary direction: 

F(v) = i [ dri- T{irf) ■ e^ z ~ D ) - 
Jo 

_ ie Mz - D) r dTj . r(cJ + i7l) . e -r)(z - D) 5 
Jo 

which, to first order in l/(z — D) is equal to 

^) = -^[3"(0)-^-%(4 
z — D 

But T(0) = 0, hence 

iuj(z - D) 

F(cj) = — ——T(u). 

i\z — D) 

Substituting this into 140(1 we obtain 

n ; i(z-D) J oj-n + iT/2 ' 

1 ,-i(fi - iT/2)h f°° , . T(w)e-Mt -* + D-h) 



V 1 i(z-D) J us-Q + iT/2 

/,(*) = --i— e -^-^/2). rj^e^m 
v 7 i(z~D) J uj-Q + iT/2 ' 

h(z) = _L_ e -*(« - *T/2)*i / dw i_i^l£ . 

v y i(z-D) J uj-Q + iT/2 

The direction of deformation of the contours in these integrals remain the same 

as it was in ff 3CJP since the factor e~ ~~ ^ in T(w) leaves the behaviour of 
the integrand at infinity unchanged. We may, therefore, write for the first two 
line of gOJl 

Ij(z) = T(Cl - iT/2) ■ I°(z - D) 
and an analogous expression with T* for the remaining lines. Hence 



|Q 3e -T(f2^r/2)| 2 a _ Ttl 
4tt(z- Df 



w{hM) = '^ e ;_r: z ■ e_i li ■ *(*2 (42) 



Since this is a sharp distribution it can be interpreted in terms of a tunneling 
time equal to zero: The velocity of tunneling is infinitely large. In the context 
of the present work this behaviour is in no conflict with the requirement that 
no information be transmittable faster than light. The reason is that the WW- 
approximation limits the validity of our calculation to ti,t% > z, i.e. to the 
inside of the light cone L of the source state preparation event, which was the 
last occasion when the experimentalist had access to the source (see Fig. 2). 
When the lifetime r of the source atom is much larger than z/c (which may 
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be of the order of several nanoseconds) the detection events fall predominantly 
within L, permitting thereby no faster than light information transfer. It is this 
region which is covered by our calculation. An improved treatment valid for 
t < z/c too (i.e. in both the vicinity of L and outside it) would certainly be of 
great interest 2 . 



Source 
world line i ^ 



iCt, Optical ct 2A Detector 
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event 



barrier 



w / 



/ 



world line 

Detection 
M event 



W: The photon "path" 

L: The light cone of the 
state preparation 
event 



State preparation * t» 
event 



Figure 2. 

For a real barrier — or even for our model barrier in a better approximation 
— the tunneling velocity will probably have a finite value which is greater than 
the velocity of light in vacuo. From the point of view of relativity theory, 
however, the point of demarcation is at c. Hence, the conlusions drawn from 
Figure 2 remain practically unchanged for any tunneling velocity larger than c. 



8 Final remarks 

Tunneling time measurements are of two very different kinds: Stationary mea- 
surements in which no moments of time are identified at all and experiments 
in which moments of time of certain real events are determined by using some 
kind of clocks. 

The Berkeley-experiment discussed in Sec.l is an ingenious example of the 
first type. No experiment of the second type has so far been performed since 
it would require precise measurement of extremely small time intervals. Even 
the purely theoretical analysis of this latter kind of experiments presents a 
challenge. The present work is an attempt to predict the result of such an 
experiment. Though the calculation performed is based on a version of the 
counting rate formulas widely used in quantum optics it cannot be considered 
completely satisfactory 11,. Quantum theory provides unambiguous rules for 

2 The limit r > would be of special importance since it is closely related to the situation 

when the photon is released at the moment of pressing a " release button" by the experimen- 
talist in any freely chosen instant of time. When r < z/c the precise nature of the state 
preparation event requires also closer examination. 
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the calculation of the probability p(t\,t2) from which the correlation function 
w(t\,t2) is obtained by differentiation. Though the rules of quantum theory 
ensure the positivity of p(ti, £2) they don't render it a nondecreasing function of 
its arguments and so the procedure may end with a negative probability density. 

The origin of this "positivity problem" may be traced back to the replace- 
ment of the spontaneous state reduction of the detectors — a process which 
falls outside the scope of the Schrodinger equation — by a "naive reduction 
hypothesis" , consisting in the identification of the statement " The detector has 
clicked" with the statement "The detector atom is in one of its excited states". 
This replacement, however, may be accepted only in the limit of weak coupling 
between the field and the detector when the latter probability is always a nonde- 
creasing function of time. In the general case the rules of quantum theory do not 
exclude the possibility that this probability decreases in some intervals of time, 
while the very notion of the "detector" is irreconcilable with such behaviour: 
For a detector the probability of being excited must never decrease. 

A possible solution of the positivity problem would be to take into account 
in the dynamics of the detector atom the influence of the equipment which it 
is built into. This would result in introducing some element of irreversibility 
into the detector's behaviour which might lead to a never decreasing excitation 
probability. Theories with spontaneous reduction (^H], might be also of 
significance in this respect. Since to first order in the detector-field interaction 
no positivity problem arises, it may, perhaps, be reasonably expected that in 
the weak coupling limit the future detector theory will be essentially reduced to 
our atomic detectors treated in first order perturbation theory on the basis of 
the naive reduction hypothesis. 

We may hence conclude that an experiment of the second kind might well 
contradict the theory in its present state even if the calculations themselves are 
irreproachable, reflecting thereby our insufficient knowledge of quantum physics 
and, perhaps, suggesting the direction toward its completion. 

The situation with the experiments of the first kind is quite the opposite. 
They belong to the domain of phenomena where the applicability of quantum 
theory has already been abundantly demonstrated. Therefore, their outcome 
can in principle be calculated in advance and the corresponding time parame- 
ter be deduced from this calculation: no contradiction with known principles 
is expected. In particular, superluminal tunneling under such stationary cir- 
cumstances never contradicts special relativity since the tunneling process is 
not accompanied by flow of information, referring to moments of time. In an 
experiment with clocks, on the other hand, superluminal tunneling would in 
general contradict relativity theory. Since the theoretical analysis of this exper- 
iment performed in the present work does not exclude completely the possibility 
of superluminal information transfer (see the end of the previous section) the 
situation deserves careful consideration. 

Acknowledgement: The author is deeply indebted to Gabor Hrasko with- 
out whose inspiring curiosity this research would not have been pursued. 
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A Derivation of the counting rate formula 

In this Appendix the derivation of the counting rate formula ©, J3J based on 
the lecture [4] is outlined. 

Consider the system, consisting of the source atom, the photon field and 
the atomic photon detector. The detector which signalizes the moment of the 
photon emission need not be considered explicitely. The total Hamiltonian 
of this system is H = H s + Hj + H s f + H c i + Hdf where H s , Hj, Hd are 
the Hamiltonians of the source, the field and the detector respectively while 
H s f 7 Hdf are the corresponding interactions. 

In order to incorporate into the calculation the irreversible nature of the 
observations of the source and the photon detector at the moments t\ and ti 
we assume that at these moments the corresponding interactions H s f and Hdf 
are switched off. This assumption will be referred to as the "irreversibility 
hypothesis" . 

The ground state and the excited states of the detector will be labelled by 
7 and e. The initial state of the system is 



Let us work in the Heisenberg picture (labelled by the superscript h) in which the 
dynamical quantities are driven by H.. After the moment of the first observation 
the state of the system becomes 



Here P g and "P e are projectors on the ground state of the source and the excited 
state e of the detector. 

After the second observation the state becomes 



where T denotes time-ordering. 

After having performed the observations the source is in its ground state, 
the detector is in one of its excited states and no photon is present. The state 
|ti, i2, e) is, therefore, identical to \g, vac, e) except that its norm is smaller than 
unity. The probability p(t\, t 2 ) introduced in Sec. 2 is equal to the square of this 
norm summed over e: 



|0) = |0) ® | 7 ) ee \e,vac,j). 




\h,h,e) 



P £ ft (i 2 )P g ?1 (i 1 )|0) if ii < t 2 , 
Pi?(h)V*(t2)\0) if i 2 < t x , 



wich can also be written as 



\t 1: t 2 ,e) =T(V?{t 2 )P*(t 1 ))\0), 




(43) 



e 
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Our aim now is to eliminate from this formula the explicite reference to the 
photon detector (except its spectral sensitivity). 

Introduce the interaction picture labelled by i by means of the unitary op- 
erator 

V{t,0) = e in o t -e- iHH , 

where s indicates Schrodinger-picture and Hq = H s — H^. In this picture the 
development of the states is governed by H % df (t) = V(t,0)H% f (t)V + (t,0): 

iV(t,0) = Hi f (t)V(t,0), 
the solution of which to first order in the detector-field interaction is 



V(t,0) = l-i / dr-HUr). 



(44) 



Since V\ = i[H 6 ,T z e } = we have V\{t) = V\ and hence 



T(V^t 2 )P^ tl )) 



V+(t2,0)VZV(t2,h)P<(it 1 )V(t 1 ,0) if t a << 2 , 
V+(t 1 ,0)Pj ! (t 1 )V{t 1 ,t 2 )T?V(h,0) if t 2 < t x . 



These expressions are to be calculated to first order, using l|44ll . 
The first line (ti < t 2 ) gives 

V+(ta,0)7>'i*(ti)- 

-iv+(t 2 ,o)v s e \( 2 dt- m f (t)p;( tl ) + / 1 dt- PlihW^t) 

Jti Jo 

= V + (t 2 ,0)V s e P;(t 1 )-iV + (t 2 ,0)V s e r dt-T(H* df (t)P l g ( tl )). 

Jo 

For the second line (t 2 < ti) we have analogously 



dt-PUt^HUt)^ 



dt-Pi(h)V s e HUt) 



In the first integral of the last line t > t 2 . Therefore, by the irreversibility 
hypothesis we have Hdf(t) = in it. Moreover, in the remaining term T>1 can 
be brought in front of the integral since P l g {ti) does not depend on H^. We 
have, therefore, 

T(rHh)P£(ti)) = C(h,h) - iV+(ti,0)V s e f 2 dt ■ T{W df (t)P l g ( tl )), (45) 

Jo 

where 

v 11 \v+(t 1 ,o)p^t 1 )V! ift 2 <h. 



18 



When l)45|l is substituted into 143|) the term Cfa, ti) gives no contribution since 
P e s |0) = 0: 



I2 



{g,vac,e\V + {t u Q)Vl / dt ■T(W df (t)P l g (t 1 ))\0) 



Now, to first order in Hm the operator V + must be replaced by unity and since 
(g, vac, e\V* = (<?, vac, e|, we have 



P 



(*i,*2)=y; (9,vac,e\ f dt-TiH^Pfr^lO) 
Jo 



Assume now that H^f — q ■ if where q acts in the Hilbert-space of the 
photon detector. For an arbitrary dynamical quantity O the z-picture and the 
Schrodinger-picture are connected by the relation 



(46) 



The Hamiltonian Ho is the sum of Hd and the Hamiltonian of the source-field 
system (H s + Hf + H s f) which commute with each other. Hence for q equation 
ffify reduces to 

q\t) = e^qe-^d*, 

and for such operators as ip and P g which are independent of the photon detector 
it gives O l = O h where the superscript H refers to the Heisenberg-picture 
introduced in Sec. 2. 

Assuming, that H^\j) — and ifj|e) = w e |e) we have 



(ti J fe)=^|(eM 7 )| 2 - (9,vac\ / dt-e w ' t T(<p H (t,P)lf(t 1 ))\0,vac) 



(47) 



Since the spectral sensitivity is given by the relation 
c^)=5>( W - We )|< e M7>| 2 , 

1471) becomes identical to the working formulae given in Sec. 2. 



B Derivation of the formula ( 1241) 

Consider the solution v + (uo z \z) of the equation JJJJ given in H10J1 . Since the 
equation is a real linear one, the combination 



v* + (ui z \z) - R*v + (u z \z) 



if z < 



a. 



. 1 - \R\ p -iw t (z-t) 
T * e iuj z (a - b) I e 

c -imjz - b) R*T c 2iu)Jb - a) c imjz - b) jf z > 5 
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is also a solution which contains an incoming wave from the right. Comparing 
this solution with (|llf> we have 

= 1-\R\ 2 

R i _ R * T c 2iuj 7 jb - a) 

The first of these equations combined with the conservation of the probability 
\T\ 2 + \R,\ 2 = 1 gives T' — T while the second one can be rewritten in the form 

TR * e -2iuj z a + T * R , e -2iu z b = Q ^ 

which is used in Sec. 5. 
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